Results from a Low-Energy Analysis of the CDMS II Germanium Data 
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We report results from a reanalysis of data from the Cryogenic Dark Matter Search (CDMS II) ex- 
periment at the Soudan Underground Laboratory. Data taken between October 2006 and September 
2008 using eight germanium detectors are reanalyzed with a lowered, 2 keV recoil-energy thresh- 
old, to give increased sensitivity to interactions from Weakly Interacting Massive Particles (WIMPs) 
with masses below ~10 GeV/c 2 . This analysis provides stronger constraints than previous CDMS II 
results for WIMP masses below 9 GeV/c 2 and excludes parameter space associated with possible 
low-mass WIMP signals from the DAMA/LIBRA and CoGeNT experiments. 



PACS numbers: 14.80. Ly, 95.35. +d, 95.30.Cq, 95.30.-k, 85.25.Oj, 29.40.Wk 



A convergence of astrophysical observations indicates 
that ^80% of the matter in the universe is in the form 
of non-baryonic, non-luminous dark matter fl. Weakly 
Interacting Massive Particles (WIMPs) Q, with masses 
from a few GeV/c 2 to a few TeV/c 2 , form a well- 
motivated class of candidates for this dark matter 0, Q . 
If WIMPs account for the dark matter, they may be de- 
tectable through their elastic scattering with nuclei in 
terrestrial detectors 

Although many models of physics beyond the Stan- 
dard Model provide WIMP candidates, supersymmet- 
ric (SUSY) models where the lightest superpartner is a 
cosmologically stable WIMP are among the most pop- 
ular 0, IH ■ In the Minimal Supersymmetric Standard 
Model (MSSM), WIMPs with masses <40 GeV/c 2 are 



generally disfavored by accelerator constraints (e.g., [j|). 
Interest in lower-mass WIMPs has been renewed by re- 
cent results from the DAMA/LIBRA @ and CoGeNT Q 
experiments, which have been interpreted in terms of 
elastic scatters from a WIMP with mass ~10 GeV/c 2 and 
cross-section ~10~ 40 cm 2 [1,0]. Although it is difficult 
to accommodate a WIMP with these properties in the 
MSSM 1 1 01] . alternate models avoid existing constraints 
(e.g.,Efr 

The CDMS II experiment attempts to identify nuclear 
recoils from WIMPs in an array of particle detectors 
by measuring both the ionization and non-equilibrium 
phonons created by each particle interaction. Back- 
grounds can be rejected on an event-by-event basis since 
they primarily scatter from electrons in the detector, de- 
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positing significantly more ionization than a nuclear re- 
coil of the same energy. Previous analyses of CDMS II 
data [lH imposed a recoil-energy threshold of 10 keV 
to maintain sufficient rejection of electron recoils that 
only ^0.5 background events would be expected in the 
signal region. At lower energies, the discrimination be- 
tween nuclear and electron recoils degrades, leading to 
higher expected backgrounds. Since WIMPs with masses 
<10 GeV/c 2 primarily produce <10 keV recoils, this 
analysis lowers the recoil-energy threshold to 2 keV, com- 
parable to the hardware trigger threshold. This lower en- 
ergy threshold increases sensitivity to low-mass WIMPs 
at the cost of significant acceptance of backgrounds. 

The data analyzed here were collected using all 30 Z- 
sensitive Ionization and Phonon (ZIP) detectors installed 
at the Soudan Underground Laboratory [13, EH ■ The de- 
tector array consisted of 19 Ge (~230 g each) and 11 Si 
(~105 g each) detectors, each a disk ~10 mm thick and 
76 mm in diameter. Each detector was instrumented 
with four phonon sensors on one face and two concentric 
charge electrodes on the opposite face. A small electric 
field (3-4 V/cm) was applied across the detectors to ex- 
tract charge carriers created by particle interactions. The 
detectors were arranged in five "towers," and are identi- 
fied by their tower number (T1-T5) and by their ordering 
within the tower (Z1-Z6). The entire array was cooled to 
<50 mK and surrounded by passive lead and polyethy- 
lene shielding. An outer plastic scintillator veto was 
used to identify showers containing cosmogenic muons 
which were not shielded by the rock overburden above 
the Soudan laboratory (2090 meters water equivalent). 

The data were taken during six data runs from Oc- 
tober 2006 to September 2008 [Hj]. Only the eight Ge 
detectors with the lowest trigger thresholds were used 
to identify WIMP candidate events since they have the 
best expected sensitivity to WIMPs with masses from 
5-10 GeV/c 2 . All 30 detectors were used to veto events 
that deposited energy in multiple detectors. 

Each detector was monitored throughout the data runs 
and periods of abnormal detector performance were re- 
moved 12j. Data taken within 20 days following ex- 



posure of the detectors to a neutron calibration source 
were removed to reduce low-energy electron-recoil back- 
grounds due to activation of the detectors. The data 
were randomly divided into two subsets before defining 
selection criteria at low energy. One subset, consisting of 
one quarter of the data (the "open" data) , was reserved 
to study backgrounds at low energy and was not used to 
calculate exclusion limits. The remaining subset totaled 
241 kg days raw exposure, after removing the bad data 
periods described above. 

The detector response to electron and nuclear recoils 
was calibrated by regular exposures of the detectors to 7- 
ray ( 133 Ba) and neutron ( 252 Cf) sources. The ionization 
energy scale was initially calibrated using the 356 keV 
line from the 133 Ba source. The phonon energy was then 



calibrated by normalizing the phonon-based recoil energy 
for electron recoils to their mean ionization energy. In 
contrast to previous analyses, a position-dependent cal- 
ibration was not applied since position-dependent varia- 
tions in the reconstructed phonon energies are less sig- 
nificant than noise fluctuations at low energies. Using 
the observed positions of the 1.3 keV and 10.4 keV acti- 
vation lines produced from exposure of the Ge detectors 
to the 252 Cf source, a small rescaling (~4%) was applied 
to ensure that the recoil energy scale for electron recoils 
was not underestimated at the 90% confidence level. For 
nuclear recoils, the recoil energy was reconstructed from 
the measured phonon energy alone by subtracting the 
Neganov-Luke phonon contribution [l4[ corresponding to 
the mean ionization measured for nuclear recoils of the 
same phonon energy from the 252 Cf source. The ratio 
of ionization to recoil energy ( "ionization yield" ) for nu- 
clear recoils was measured down to ^4 keV, below which 
a power-law extrapolation was used. 

Candidate events were required to pass basic recon- 
struction quality cuts similar to the criteria used in pre- 
vious analyses of these data [12J. Due to the negligible 
probability of a WIMP interacting more than once in 
the apparatus, candidates were required to have energies 
consistent with noise in all but one detector and have no 
coincident activity in the plastic scintillator veto. They 
were further required to have ionization signals consistent 
with noise in the outer charge electrode. The ionization 
energy was required to be within (+1.25, — 0.5)ct of the 
mean ionization energy for nuclear recoils measured from 
calibration data, which defines the "nuclear-recoil band." 
This asymmetric band, which has been tig htened rela- 
tive to previous low-energy analyses IjJ was chosen 
based on calibration data and the observed low-energy 
backgrounds in the open data in order to maximize sensi- 
tivity to nuclear recoils while limiting leakage from elec- 
tron recoils and zero-charge events. The recoil-energy 
range considered for this analysis was 2-100 keV. 

The hardware trigger efficiency was determined using 
events for which at least one other detector triggered, 
which provide an unbiased selection of events near thresh- 
old. The data are well described by an error function, 
with the mean trigger threshold varying from 1.5-2.5 keV 
for the eight Ge detectors. Based on the selection criteria 
above, the signal acceptance was measured using nuclear 
recoils from the 252 Cf calibration data. We calculated 
the nuclear-recoil band acceptance conservatively by as- 
suming all events with ionization energy < 3<r above the 
mean of the distribution were nuclear recoils. In par- 
ticular, the zero-charge events described below were in- 
cluded, although their rate in the 252 Cf calibration data 
is negligible. The livetime-weighted average of the indi- 
vidual detector selection efficiencies is shown in the inset 
of Fig. [I] with the largest loss of efficiency coming from 
the requirement on ionization energy. 

The energy spectrum for the candidate events passing 
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FIG. 1. (color online). Comparison of the energy spectra 
for the candidate events and background estimates, co-added 
over the 8 detectors used in this analysis. The observed event 
rate (error bars) agrees well with the electron-recoil back- 
ground estimate (solid), which is a sum of the contributions 
from zero-charge events (dashed), surface events (+), bulk 
events (dash-dotted), and the 1.3 keV line (dotted). The se- 
lection efficiencies have been applied to the background es- 
timates for direct comparison with the observed rate, which 
does not include a correction for the nuclear-recoil acceptance. 
The inset shows the measured nuclear-recoil acceptance effi- 
ciency, averaged over all detectors. 

all selection cuts is shown in Fig. [TJ Although the shape 
of the observed spectrum is consistent with a WIMP sig- 
nal, we expect that a significant number of the candidates 
are due to unrejected electron recoils. Figure [2] shows 
the distribution of candidates in the ionization-yield ver- 
sus recoil-energy plane for T1Z5. Several populations of 
events which can leak into the signal region at low energy 
are apparent. For each population described below, we 
measure the rate and energy spectrum in sidebands where 
the contribution from low-mass WIMPs would be negligi- 
ble, and extrapolate the observed spectrum to lower ener- 
gies to estimate the leakage. The systematic errors intro- 
duced by these extrapolations are potentially large and 
are difficult to quantify. However, as shown in Fig. [T] and 
discussed below, these simple extrapolations can plausi- 
bly explain all the observed candidates. 

Events with ionization energies consistent with noise 
are seen below the nuclear-recoil band. Most or all 
of these "zero-charge" events arise from electron recoils 
near the edge of the detector, where the charge carri- 
ers can be completely collected on the cylindrical wall 
rather than on the readout electrodes. At recoil energies 
>10 keV, these events can be rejected using a phonon- 
based fiducial- volume cut. At lower energies, reconstruc- 
tion of the event radius using phonon information is un- 
reliable. To maintain acceptance of low-energy nuclear 
recoils, some zero-charge events are not rejected at ener- 
gies <5 keV where the ionization signal for nuclear recoils 
becomes comparable to noise. By extrapolating the expo- 
nential spectrum observed for zero-charge events above 
5 keV, we estimate that they contribute ~50% of the 
candidate events. 
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FIG. 2. (color online). Events in the ionization-yield 
versus recoil-energy plane for T1Z5. Events within the 
(+1.25, — 0.5)cr nuclear-recoil band (solid) are WIMP candi- 
dates (large dots). Events outside these bands (small, dark 
dots) pass all selection criteria except the ionization-energy 
requirement. The widths of the band edges denote variations 
between data runs. Events from the 252 Cf calibration data 
are also shown (small, light dots). The recoil-energy scale as- 
sumes the ionization signal is consistent with a nuclear recoil, 
causing electron recoils to be shifted to higher recoil energies 
and lower yields. 

A second source of misidentified electron recoils comes 
from events interacting near the detector surfaces, where 
ionization collection may be incomplete. These events 
are primarily concentrated just above the nuclear-recoil 
band, with an increased fraction leaking into the sig- 
nal region at low energies. For recoil energies >10 keV, 
nearly all such surface events can be rejected [l2j be- 
cause they have faster-rising phonon pulses than nuclear 
recoils in the bulk of the detector. This analysis does 
not use phonon timing to reject these events since the 
signal-to-noise is too low for this method to be effective 
for recoil energies <5 keV. Extrapolating the exponen- 
tial spectrum of surface events identified above 10 keV 
implies that ~15% of the candidates are surface electron 
recoils. 

At recoil energies <5 keV, the primary ionization- 
based discrimination breaks down as the ionization sig- 
nal becomes comparable to noise even for electron recoils 
with fully collected charge. Extrapolating the roughly 
constant electron-recoil spectrum observed above 5 keV 
indicates that ~10% of the observed candidates arise 
from leakage of this background into the signal region. 
Just above threshold, there is an additional contribution 
to the constant electron-recoil spectrum from the 1.3 keV 
line, which leaks above the 2 keV analysis threshold since 
our recoil-energy estimate assumes the ionization signal 
is consistent with a nuclear recoil. The measured in- 
tensity of this line at ionization yields above the signal 
region indicates that the 1.3 keV line accounts for ~10% 
of the observed candidates. T1Z5 has less expected leak- 
age from these fully-collected electron-recoil backgrounds 
than the average detector since it has the best ionization 
resolution. 
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FIG. 3. (color online). Top: comparison of the spin- 
independent (SI) exclusion limits from these data (solid) to 
previous results in the same mass range (all at 90% C.L.). 
Limits from a low-threshold analysis of the CDMS shallow- 
site data El (dashed), CDMS II Ge results with a 10 keV 
threshold [12] (dash-dotted), recalculated for lower WIMP 
masses, and XENON100 with constant (+) or decreasing (□) 
scintillation-efficiency extrapolations at low energy [17J are 
also shown. The filled regions indicate possible signal regions 
from DAMA/LIBRA (dark), CoGeNT (light) and 
a combined fit to the DAMA/LIBRA and CoGeNT data @] 
(hatched). Bottom: comparison of the WIMP-neutron spin- 
dependent (SD) exclusion limits from these data (solid), 
CDMS II Ge results with a 10 keV threshold (dash-dotted), 
XENON10 [Hi (A), and CRESST l?] (O). The filled re- 
gion denotes the 99.7% C.L. DAMA/LIBRA allowed region 
for neutron-only scattering [20|| . An escape velocity of 544 
km/s was used for the CDMS and XENON100 exclusion lim- 
its, whereas the other results assume an escape velocity from 
600-650 km/s. Using the same halo parameters as assumed 
for the allowed regions would lead to slightly stronger limits 
(dotted). 



These estimates indicate that we can claim no evi- 
dence for a WIMP signal. However, since the back- 
ground model involves sufficient extrapolation that sys- 
tematic errors are difficult to quantify, we do not sub- 
tract backgrounds but instead set upper limits on the 
allowed WIMP-nucleon scattering cross section by con- 
servatively assuming all observed events could be from 
WIMPs. Limits are calculated using the high statistics 
version of Yellin's optimum interval method [2l[ . Data 
from multiple detectors are concatenated as described 
in [l5| . This method allows the choice of the most con- 
straining energy interval on the lowest background de- 
tector while applying the appropriate statistical penalty 



for the freedom to choose this interval. The method and 
the ordering of detectors by position within the tower 
were specified with no knowledge of the WIMP candi- 
dates to avoid bias. For WIMP masses from 5-8 GeV/c 2 , 
the most constraining interval contains events only from 
T1Z5 and has no dependence on the detector ordering 
used. The standard halo model described in [22| is used, 
with specific parameters given in 15 , 23| . The candidate 
event energies and selection efficiencies for each detector 
are given in [23j . 

The limits do not depend strongly on the extrapola- 
tion of the ionization yield used at low energies since the 
Neganov-Luke phonon contribution is small for recoil en- 
ergies below 4 keV. Conservatively assuming 25% lower 
ionization yield near threshold would lead to only ^5% 
weaker limits in the 5-10 GeV/c 2 mass range. 

Figure [3] (upper panel) shows the resulting 90% up- 
per confidence limit on the spin-independent WIMP- 
nucleon scattering cross section. This analysis provides 
stronger limits than previous CDMS II Ge results for 
WIMP masses below ~9 GeV/c 2 , and excludes param- 
eter space previously excluded only by the XENON10 
and XENON100 experiments for a constant extrapola- 
tion of the liquid xenon scintillation response for nuclear 
recoils below 5 keV 17|, |24[, |25[ . Our analysis provides 
stronger constraints than XENON10 and XENON100 be- 
low ^7 GeV/c 2 under conservative assumptions for the 
scintillation response [H, [T3, HH . 

Spin-dependent limits on the WIMP-neutron cross sec- 
tion are shown in Fig. [3] (lower panel), using the form fac- 



tor from [27J. XENON10 constraints, calculated assum- 



ing a constant extrapolation of the scintillation response 
at low energy lj| 26[ , are stronger than these results for 
WIMP masses above ~7 GeV/c 2 . 

These results exclude interpretations of the 
DAMA/LIBRA annual modulation signal in terms 
of spin-independent elastic scattering of low-mass 
WIMPs (e.g., d [ll]). We ignore the effect of ion 
channeling on the DAMA/LIBRA allowed regions since 
recent analyses indicate channeling should be negligi- 
ble [ill [28[ . These results are also incompatible with a 
low-mass WIMP explanation for the low-energy events 
seen in CoGeNT @,tl]. 

The CDMS collaboration gratefully acknowledges the 
contributions of numerous engineers and technicians; we 
would like to especially thank Jim Beaty, Bruce Hines, 
Larry Novak, Richard Schmitt and Astrid Tomada. In 
addition, we gratefully acknowledge assistance from the 
staff of the Soudan Underground Laboratory and the 
Minnesota Department of Natural Resources. This work 
is supported in part by the National Science Founda- 
tion (Grant Nos. AST-9978911, PHY-0542066, PHY- 
0503729, PHY-0503629, PHY-0503641, PHY-0504224, 
PHY-0705052, PHY-0801708, PHY-0801712, PHY- 
0802575, PHY-0847342, and PHY-0855525), by the De- 
partment of Energy (Contracts DE-AC03-76SF00098, 



5 



DE-FG02-91ER40688, DE-FG02-92ER40701, DE-FG03- 
90ER40569, and DE-FG03-91ER40618), by the Swiss 
National Foundation (SNF Grant No. 20-118119), and 
by NSERC Canada (Grant SAPIN 341314-07). 

APPENDIX 

The following appendix includes additional details re- 
garding the calibration of the energy scale and expected 
backgrounds to clarify several of the points raised in [2i| . 

The phonon-based recoil energy scale for this analy- 
sis was calibrated using the position of the 1.298 keV 
and 10.367 keV activation lines which conveniently lie in 
the recoil energy range of interest for low mass WIMPs. 
These lines provide a robust determination of the en- 
ergy scale for electron recoils near threshold. We use 
the most conservative values for the energy scale which 
are consistent with the position of these activation lines 
at the 90% CL for each detector. Figure |4] shows the 
measured positions of the 1.3 keV and 10.4 keV acti- 
vation lines for T1Z5, after calibrating the energy scale 
to ensure that the energy is not underestimated for these 
lines. The recoil energy units are keVee, or keV "electron- 
equivalent" , which indicates that the total phonon signal 
has been corrected for the contribution from Neganov- 
Luke phonons assuming that the ionization produced was 
consistent with an electron recoil. 

Using the phonon energy scale calibrated with the 
electron-recoil activation lines, the recoil energy for the 
candidate nuclear recoils was calculated. This recoil en- 
ergy estimate uses only the phonon signal and does not 
use the ionization signal on an event-by-event basis due 
to the poorer signal-to-noise of the ionization pulses. The 
measured total phonon signal is corrected for the small 
contribution from Neganov-Luke phonons (~15% of the 
phonon signal for low-energy nuclear recoils) by subtract- 
ing the drift heat corresponding to the mean ionization 
energy measured for nuclear recoils in the 252 Cf calibra- 
tion data. Figure [5] shows the ionization energy versus 
recoil energy for the WIMP search events and 252 Cf cal- 
ibration data, using the recoil energy determined by as- 
suming the Neganov-Luke phonon contribution is consis- 
tent with a nuclear recoil. 

To properly correct for the Neganov-Luke phonons, we 
must determine the ionization yield for nuclear recoils as 
a function of recoil energy. Above ^4 keV, depending on 
detector, the peak of the nuclear recoil ionization distri- 
bution lies > 2<j above noise, and the ionization yield was 
found by fitting this distribution to a Gaussian in bins 
of recoil energy. A power law was then fit to these yield 
data from ^4-20 keV and extrapolated to determine the 
ionization yield at lower energies, as shown in Fig. |6] (left) 
for T1Z5. 

As a cross-check on this extrapolation, we can extend 
the fits to lower energies by fitting only the portion of the 
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FIG. 4. Efficiency-corrected phonon recoil-energy spectrum 
for electron recoils in T1Z5. The solid lines show fits to the 
location of the activation lines, which give mean values of 
1.333±0.025 keVee and 10.391±0.022 keVee. The resolution 
of the 1.3 keV line is ~100 eVee, consistent with the expected 
resolution from noise. The 10.4 keV line is broadened by posi- 
tion dependence of the phonon signal for which no correction 
has been applied. Cosmogenic 65 Zn has decayed away suf- 
ficiently that its contribution to the low energy tail of the 
activation peaks is negligible [3C| . The relative intensity of 
the lines measured from the fits is 0.145±0.030, consistent 
with expectations. 
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FIG. 5. Events in the ionization energy vs. recoil energy plane 
for T1Z5. Events from the 252 Cf calibration data (small, gray 
dots) and WIMP search data (large, black dots) are shown. 
The recoil energy scale is given by the total phonon energy mi- 
nus the Neganov-Luke phonon contribution corresponding to 
the mean ionization for nuclear recoils. This scale, in units of 
keVnr, gives the correct recoil energy only for nuclear recoils, 
while electron recoils appear at higher recoil energy due to 
the larger contribution of Neganov-Luke phonons. The solid 
lines show fits to the mean of the nuclear recoils (green) and 
electron recoils (blue). The red dotted curves indicate lines 
of constant recoil energy when the ionization signal is used to 
determine the Neganov-Luke phonon contribution. 
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FIG. 6. (left) Measurement of the mean ionization yield for nuclear recoils in T1Z5. Both 3-parameter fits to the yield 
distribution (black, squares) and 1-parameter fits where the width and amplitude are constrained (blue, circles) are shown. 
The black solid line shows the power-law fit to the measured yields above 4 keV, which was extrapolated to lower energies. 
This fit was based on data which improperly included events outside the fiducial volume, leading to slightly lower yields than 
shown by the data points from ~7-15 keV. This shift has a negligible effect on the energy scale. Several other models for the 
ionization yield at low energy are also shown: Lindhard prediction [22l . I3l1 ] (red, dashed), steeper yield extrapolation below 
4 keV (magenta, dotted), and Lindhard, with fc=0.1 (green, dot-dashed), (right) Effect on the exclusion limits from this analysis 
for the ionization yield models shown on the left. The limits are only weakly dependent on the extrapolation of the ionization 
yield at low energies since the Neganov-Luke phonon contribution is small for low-energy nuclear recoils. The 90% CL allowed 
regions for CoGeNT (blue), DAMA/LIBRA (gray), and a combined fit to CoGeNT and DAMA/LIBRA (hatched) from |j| 
are shown. The region where models including light WIMPs were found to provide a good fit to CoGeNT in 0] is also shown 
(orange). The limits were calculated using i>o=220 km/s and u esc =544 km/s, with detailed parameters given in [2^ ]. 



distribution which lies > 2a above noise. To ensure that 
these fits are well-behaved at energies where the peak of 
the distribution is not contained in the fitting window, we 
constrain the width of the Gaussian using the measured 
ionization resolution as a function of energy from acti- 
vation lines at 1.3, 10.4 and 66.7 keV. The amplitude of 
the Gaussian is fixed based on Monte Carlo simulations of 
the 252 Cf neutron calibration recoil-energy spectrum. We 
then perform a 1-parameter fit to determine the mean of 
the distribution. Both fitting methods give similar yields 
above ~4 keV, while the fixed width and amplitude fits 
can be extended down to the threshold of 2 keV. The 
yields determined by the 1-parameter fits agree well with 
the power-law extrapolation. Details of the measurement 
of the ionization yield will be presented in an upcoming 
publication [32| . 

Figure H] (right), shows that the energy scale and lim- 
its are only weakly dependent on the extrapolation of 
the ionization yield at low energy since the Neganov- 
Luke phonons contribute only a small fraction of the total 
phonon signal. Conservatively assuming ionization yields 
which are ~25% lower than those suggested by our data, 
and by previous measurements by other groups in this 
energy range (see e.g. [|| and references therein), would 
not change the conclusions of this analysis. 

Although the correction for the Neganov-Luke phonons 
does not contribute significant uncertainty to the recoil 



energy scale, larger errors are possible if the phonon 
collection for nuclear recoils differs from that for elec- 
tron recoils. The absolute nuclear recoil energy scale is 
constrained by the measured ionization yields shown in 
Fig. [Bl which are ~15% lower than previous measure- 
ments in the energy range of interest. Assuming that the 
ionization produced is consistent with previous measure- 
ments of the quenching factor, the measured yields imply 
an overestimate of the recoil energy scale near thresh- 
old. To be conservative, we do not apply a correction to 
the recoil energy scale based on the measured yields, but 
such a correction would improve our constraints in the 
5-10 GeV/c 2 mass range. 

We calculated limits on the WIMP-nucleon cross sec- 
tion by conservatively assuming all the candidate events 
could arise from WIMP-induced nuclear recoils. How- 
ever, we expect significant backgrounds at low energy 
which can mimic a WIMP signal. Although estimates of 
these backgrounds do not affect the limits on the WIMP- 
nucleon cross section derived, they indicate that we can- 
not claim evidence for a WIMP signal since all observed 
candidates can plausibly be accounted for by expected 
backgrounds. 

In particular, zero-charge events are expected to be 
the limiting background in the energy range of inter- 
est for WIMPs with masses ~7 GeV/c 2 . These events 
are consistent with electron recoils at very high radius 
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FIG. 7. Zero-charge event rate vs. recoil energy. The ob- 
served rates for single- and multiple-scatter zero-charge events 
in the WIMP search data are shown, both for the total rate 
coadded over detectors and for T1Z5 alone. Below ~5 keV, 
bulk electron recoils can leak into the zero-charge band, lead- 
ing to a more rapid increase in the observed spectrum. The 
rate of zero-charge events in the 133 Ba calibration data is also 
shown, after subtracting the expected contribution at low en- 
ergy from the constant bulk electron-recoil spectrum and scal- 
ing by a factor of ~5 x 10 -4 to match the zero-charge event 
rate in the WIMP search data from 5-10 keV. The solid line 
shows the exponential fit to the coadded WIMP search singles 
rate from 5-15 keV, which was extrapolated to estimate the 
zero-charge event rate near threshold. 



in the detector, where the ionization can be completely 
collected on the cylindrical walls instead of on the ion- 
ization electrodes. Zero-charge events are seen in the 
WIMP search data for both events interacting in a sin- 
gle detector ( "singles" ) and events with energy deposited 
in more than one detector ("multiples"). As shown in 
Fig. El the zero-charge multiples and singles are observed 
to occur with similar rates and approximately exponen- 
tial spectra below ~15 keV. Due to the negligible prob- 
ability of a WIMP interacting in more than one detec- 
tor, the multiple-scatter zero-charge events can arise only 
from backgrounds. We expect a corresponding contribu- 
tion to the zero-charge singles rate from the same back- 
ground sources when energy is deposited in inactive por- 
tions of the experiment and only a single detector. For 
bulk electron recoils, the rate of tagging multiple scatters 
is approximately twice the single-scatter rate. For events 
interacting on the cylindrical surfaces of the detectors, 
the multiple-scatter tagging is expected to be less effi- 
cient since the cylindrical walls face the copper detector 
housings rather than another detector. 

In addition, zero-charge events are seen in the 133 Ba 
calibration data with a similar spectrum to the WIMP 
search data. Their rate in the 133 Ba calibration data 
is significantly higher, with ~4 times more zero-charge 
events observed from 5-10 keV in an exposure which 



was ^1/500 as large. This demonstrates that electron 
recoils can produce zero-charge events with an increas- 
ing spectrum at low energy. Although the low-energy 
gamma spectrum is approximately constant, we expect 
that zero-charge events should increase in rate at low en- 
ergy since only those events which interact in a small 
dead layer near the cylindrical wall will have their charge 
trapped and pass the fiducial volume selection. Due to 
the rapidly decreasing penetration depth of low-energy 
external gammas, an increasing fraction of low-energy 
events will interact in this dead layer. However, the 
zero-charge event rate does not scale directly with the 
electron-recoil rate since only ~l/6 as many zero-charge 
events were observed in the 133 Ba calibration data as in 
the WIMP search singles and multiples from 5-10 keV, 
when normalized to the total number of bulk electron 
recoils. 

To estimate the contribution of zero-charge events at 
low energy, we extrapolated the approximately exponen- 
tial spectrum observed for the coadded WIMP search sin- 
gles from 5-15 keV to lower energies, as shown in Fig. [7] 
This extrapolation was found to be relatively insensi- 
tive to the choice of lower boundary for the fit, with a 
20% variation in the expected rate from 2-5 keV when 
the lower edge of the fitting window was varied from 5- 
8 keV. Near 5 keV, there is some overlap of the nuclear- 
recoil and zero-charge distributions. If WIMP-induced 
nuclear recoils are present, a small fraction could appear 
in the zero-charge band and bias the background esti- 
mate. Even conservatively assuming that all events in the 
(+1. 25,-0. 5)ct nuclear-recoil acceptance region are due 
to WIMP-induced nuclear recoils, we would expect only 
~15% of the observed single-scatter zero-charge events 
from 5-6 keV would be due to WIMPs. The overlap 
of the distributions decreases rapidly with energy, and 
any error introduced is small relative to possible errors 
from the extrapolation itself. Although these estimates 
indicate that we expect a significant zero-charge contri- 
bution in the region of interest for light WIMPs, we do 
not subtract this background due to possibly significant 
systematic errors introduced by the extrapolation which 
are difficult to quantify. 

A direct comparison of the CDMS and CoGeNT spec- 
tra is shown in Fig. |SJ Both the CDMS data coadded 
over all detectors and the data for T1Z5 alone are shown. 
T1Z5 has the best ionization resolution, leading to the 
best discrimination against electron-recoil backgrounds, 
and it dominates the sensitivity of this analysis in the 
5-10 GeV/c 2 mass range. Even with no background sub- 
traction, the CDMS data for T1Z5 are incompatible with 
a low-mass WIMP interpretation for the entire CoGeNT 
excess. Figure [8] also shows the expected WIMP spectra 
for the models considered in [29[ . The solid line shows the 
WIMP model described as the best fit to the CoGeNT 
data for m x —7 GeV/c 2 . As shown in Fig. [5J this model 
requires a significant fraction of the CoGeNT low-energy 
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FIG. 8. (left) Comparison of the observed rate in CoGeNT (red, squares), CDMS coadded over all detectors (blue, circles), 
and CDMS detector T1Z5 (black, triangles), which sets the strongest constraints in the 5-10 GeV/c 2 mass range. All data are 
corrected for the nuclear-recoil acceptance efficiencies. The CoGeNT data shown have the L-shell electron capture peaks and 
a constant background subtracted 29]. The CDMS ionization spectra are based on the recoil energy reconstructed from the 
phonon signal alone and have been converted to ionization energy using the ionization yields shown as the solid black curve in 
Fig El (right) Comparison of the same data versus recoil energy. The CoGeNT data have been converted to recoil energy using 
the quenching factor assumed in [8J. This quenching factor is slightly higher than the ionization yields measured by CDMS, 
causing the spectra to appear more compatible than when plotted versus ionization energy. The recoil spectra for the WIMP 
models considered in Fig. 1 of ref. [2t| are also shown. The dotted line indicates the expected spectrum for m x =7 GeV/c 2 and 
<7 SJ =1.4x 1CT 40 cm 2 , corresponding to a WIMP model from a simultaneous fit to the DAMA/LIBRA and CoGeNT data [sl,!^. 
The solid line shows the spectrum for the same WIMP mass and asj=5 x 10 -41 cm 2 , which is described as the best fit for a 
m x =7 GeV/c 2 WIMP to the CoGeNT data in @. This point lies outside the CoGeNT allowed regions from 0,11] and is not 
excluded by our analysis. The WIMP spectra assume vo=220 km/s and v eS c=544 km/s. 



excess to be due to an exponential background, and it 
is not excluded by this analysis. To make the observed 
spectra for CDMS and CoGeNT compatible would re- 
quire the majority of the CoGeNT excess to be due to 
low-energy backgrounds and the CDMS backgrounds to 
be smaller than expected. 
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